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ABSTRACT 

We investigate structures of hybrid stars, which feature quark core surrounded by 
a hadronic matter mantle, with super-strong toroidal magnetic fields in full general 
relativity. Modeling the equation of state (EOS) with a first order transition by bridg- 
ing the MIT bag model for the description of quark matter and the nuclear EOS by 
Shen et al., we numerically construct thousands of the equilibrium configurations to 
study the effects of the phase transition. It is found that the appearance of the quark 
phase can affect distributions of the magnetic fields inside the hybrid stars, making the 
maximum field strength about up to 30 % larger than for the normal neutron stars. 
Using the equilibrium configurations, we explore the possible evolutionary paths to 
the formation of hybrid stars due to the spin-down of magnetized rotating neutron 
stars. We find that the energy release by the phase transition to the hybrid stars is 
quite large (< 10 52 erg) even for super strongly magnetized compact stars. Our results 
suggest that the strong gravitational-wave emission and the sudden spin-up signature 
could be observable signals of the QCD phase transition, possibly for a source out to 
Megaparsec distances. 

Key words: dense matter - equation of state - gravitation - stars: rotation - 
stars:neutron 



1 INTRODUCTION 

A very hot issue in hadronic and nuclear physics is to 
search the phase transition from baryons to their consti- 
tutes - deconfined quarks. Heavy ion colliders such as RHIC 
(Brookhaven) and LHC (CERN) are now on line to explore 
the QCD phase diagram for the high temperature and small 
baryon density regimes, for which lattice QCD calculations 
predict a smooth crossover to the QCD phase transition (see 
IStephanovl (|2004h for review). 

Conversely for the low temperature and high baryon 
density regimes, compact stars are expected to provide a 
unique window on the phase transition at their extreme den- 
sity with super-strong magnetic field. It has been suggested 
long ago that quark matter may exist in the interio r of com- 
pact objects l|ltohlll970l : lBodmerlll97ll : IWittenll98i V Hybrid 
stars (and strange quark stars) are considered to be such ob- 
jects, which feature quark cores surrounded by a hadronic 
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matte r mantle (or quark co res only) (for reviews, e.g., Weber 
1 19991 ) ; lGlendenninel (|200ll )). Even if the relevant conditions 
could be reached in a laboratory in the near future l|Sengerl 
12004 ) . the conditions prevailing in the compact stars are dif- 
ferent from those produced in accelerators, i.e., the matter is 
long-lived, charge neutral and in /3-equilibrium with respect 
to weak interactions. It is therefore important to investigate 
the properties of such "exotic" stars, providing hints about 
the main features of matter at those extreme conditions. 

The formation of the quark cores in compact stars 
is expected to take place by a first-order phase transition 
( Glcndcn ning||l992l . 120011 ). Albeit still very uncertain (e.g., 
iHorvathl i|2007l )), such a transition would proceed by the 
conversion of initially metastable hadronic matter in the 
core into the new deconfined quark phase. The metastable 
phase could be formed as the central density of neutron stars 
exceeds a critical value, due to mass accretion, spin-down 
or cooling. Possible astrophysical cites are in the protoneu- 
tron stars during th e collapse of of supernova cores (near 
the e poch of bounce (|Takahara fc Satoll 19851: lYasutake et al.l 
I2OO7I ) or at the late postbounce phase ( Gentile et al. 19931 : 



2 N. Yasutake, K. Kiuchi, and K. Kotake 



iNakazato etafl 120081 ; ISagert et all 120081 )') and in old neu- 
tron s tars accreting from their com panions (|Benvenuto et all 
1 19941 : IChau1ll997l : iLin et alJl200d ). In whichever cases, the 
sudden nucleation of the exotic phase in the hadronic 
star will be accompanied by a core-quake and huge en- 
ergy release of the gravitational binding energy. Such en- 
ergy release has been proposed t o explain the central en- 
gines of the gamma- ray bursts (Bo mbaci fc Dattal |2000| ; 
iBerezhiani et al.ll2003l ). Possible observables of these tran- 
sient phenom enon should be glit ches, magnetar flares, and 
superbursts (|Alford et al.l l2007h. In addition, the detec- 
tion of gr a vitati o nal waves (Ilokall200ll; lYasutake et al. I l2007l ; 



ILin et al 



2006; Abdikamalov et al.l 120081) and neutrinos 
1 Nakazato et alj 120081 ; ISagert et alj 200sf ) generated at the 



moment of the phase transition, should supply us implica- 
tions to unveil the mechanism of the phase transition. 

Here it should be noted that most of these cal- 
culations/estimations concerning the phase transition 
inside compact stars are limited to a non-rotating case, 
in which the Tolman-Oppenheimer -Volkoff equat i on is 
solved to obtain t heir structures ( Haensel et ap [1 986; 
Zduniketafl Il987l; iMuto fc Tatsumil Il99d ; i Drago et al.1 



20041; I Yasutake fc Kashiwal |2009 | ) Exceptions are for 

Gourgoulhon et all (|l999l ): lYasutake et"aH (|2005h ; 
Zdunik et al.1 ( 20071 ) . in which relativistic equilibrium 
configurations of rotating strange stars, beyond the so- 
called slow rotati on approximations (see references in 
IGlendennind (1200 ll)). are constructed for e st imating the en- 
ergy r elease. In Gourgoulhon et ah! l| 19991 ); lYasutake et al.l 
l|2005l ). hadron matter is assumed to be converted fully to 
quark matter, leading to t h e form a tion of strange quark 
stars li ke in Alcock et al.l (Il986l ); Benvenuto fc Horvatbl 
l| 19891 ); lOlesen fc Madsenl (|l99ll ); iLugones et al.l l| 19941 ). 
However, it is recently pointed out that the strange quark 
stars could be ruled out by their too fast spin-down rates via 
gravi tational radiations from r-modes instability l|Madsenl 
2000). The quasi-periodic oscillations (QPOs) of strange 
quark stars could not be reconciled with the obse rv ations 
l| Watts fc Reddvl 120071 ). Moreover, iMallick et all (120091 ) 
has recently claimed that magnetars cannot convert to 
purely quark stars, but only to hybrid stars. These suggest 
us to pay attention to hybrid stars rather than stars 
ma de of purely d econfi ned quarks. In a very recent work 
by IZdunik et al.l l|2008l ). the phase transition of rotating 
hybrid stars is discussed, however, the equation of state 
(EOS) is made very idealistic, assumed to model a strong 
phase transition, seemingly less sophisticated than the 
EOSs in the recent literature cited above. In addition to 
rotation, neutron stars observed in nature are magnetized 
with the typical magnetic fi eld strength ~ 10 11 -10 13 G 
(|Lvne fc Graham-Smith] 120051 ). The field strength is often 
much larger than the canonical value as ~ 10 15 G for 
a special class of the neu tr on stars such as magnetar s 
i|Lattimer fc Prakashl 120071 ; I Woods fc Thompson! 20061) 
In a series of our recent papers dKiuchi fc Kotakef 20081 ; 
iKiuchi fc Yoshidall200Sl ; lKiuchi et al.ll2008dl) . wc have stud- 
ied equilibrium configurations of relativistic magnetized 
compact stars to apply for the understanding of their 
formations and evolutions, but with hadronic EOSs. 

Above situations motivate us to investigate the struc- 
tures of relativistic hybrid stars with magnetic fields, and 
by using them, to estimate the energy release at the phase 



transition . This study i s pose d as an extension to the 
study by lYasutake et all l|2005l ). in which the phase tran- 
sition to the rotating strange stars was investigated. Due 
to the unavailability of the method to construct a fully 
general relativistic star with arbitrarily magnetic structures 
(namely both with toroidal and poloidal fields), we here 
consider the equil i brium with purely toroidal fields as in 
IKiuchi fc Yoshidal ((2008) . It is noted th at the outcomes o f 
the recent stellar evolution calculations dHeger et al.ll2005l ) 
and the MHD(magne tohydrodynam i cs) simulations of core- 



collap se supernovae ([Kotake et al 



f ibH); ITakiwaki et al 
tessart et al.1 j|200' 
ISawai et all i|200Sft ; IKiuchi et al.1 l|2008cl ); ITakiwaki et al 



(2004 ); lObergaulinger et al.1 (|2006[) 



(2009)), suggesting much dominance of the toroidal fields, 
are not in contradict ion with the assumpti on. Following the 
scenario proposed in lYasutake et al.l(l2005l ). we consider the 
possible evolutionary tracks of a rapidly rotating and mag- 
netized neutron star to a slowing rotating hybrid star due to 
the spin-down via gravitational radiation and/or magnetic 
braking. During the evolutions, the baryon mass and the 
magnetic field strength are taken to be constant for simplic- 
ity. The energy release can be estimated from the difference 
in the mass-energies between the hadronic star and the hy- 
brid star along each sequence. By constructing thousands of 
equilibrium configurations, we hope to clarify the possible 
maximum energy release at the moment of the transition 
and discuss their astrophysical implications. 

The paper is organized as follows. The method for con- 
structing the EOS with the phase transition and the nu- 
merical scheme for the stellar equilibrium configurations are 
briefly summarized in Sec. [2] Sec. [3] is devoted to showing 
numerical results. Summary and discussion follow in Sec. [4] 
In this paper, we use geometrical units with G — c = 1. 



2 EQUATION OF STATE AND NUMERICAL 
METHOD 

2.1 Equations of State with a first-order phase 
transition 

As mentioned, we assume that the deconfinement of the 
quarks takes place at the first order phase transition. In this 
case, a mixed phase can form and it is typically described us- 
ing two separate EOSs, one for the hadronic and the other 
for the quark phase. Here the bulk Gibbs construction is 
used to bridge the two phases. In the following, we briefly 
summarize the adopted EOSs for each phase and explain 
the features of the resulting EOS with the first order phase 
transition. 

Since lattice QCD is yet to make solid predictions for 
large density regimes in compact objects, the quark mat- 
ter EOS is currently computed using phenomenological de- 
scriptions such as the MIT bag or the Nambu-Jona-Lasinio 
(NJL) models. For the quark phase, we choose the EOS 
base d on th e very si mple but wid e ly app lied MIT bag model 
(see IWeberl l| 19991 ); IGlendennind (|2001l ) for reviews). Using 
the model, one can express the energy density, e, and the 
pressure, P of strange quark matter as functions of baryon 
number density, n in the following form, 



(i) 
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P 



2 d(e/n) 



1 + x'j — arsinh a; /] 



d n 



(2) 
(3) 



where m/ is the quark mass of the flavor of / taken to be 



m d 



5 MeV and m s 



150 MeV, x; = k { P /m f 



( f) 

is a normalized Fermi wave number of kp , and B is the 
bag constant. Baryon number density of strange matter is 
n — \{n u + rid + n s ), where n; = 3 /n 2 is the num- 
ber density of / quarks. We use a si mple MIT bag mode l 
of self-bound strange quark matter (IChodos et al.l [l974h . 
neglecting the quark interactions except for the confine- 
ment effec ts described by the b ag constant (see, e.g., ref- 
erences in iBonanno et al.l (|2007l ) for sophistication of the 
model). We set that the bag constants are 200 MeV fm~ 3 
and 250 MeV mi -3 in this paper. These values seem con- 
sisten t with the implicat ions in the recent lattice QCD re- 
sults (llvanov et al.ll2005l ). The quark EOS with lower bag 
constant than 200 MeV fm -3 may be too soft to be com- 
patible with the observation of the massive neutron star like 
ulsars Ter 5 I and J (> 1.68Mr with 95 % confidence) 
Ransom et all |2005l : I Afford et all |2007| ). For the EOS of 
th e hadron phase, w e adopt the nuclear EOS developed 
by IShen et all (|l998h . which are often emplo yed in recent 
MHD simulations relevant for magnetars (see iKotake et al.l 
(2006) for a review). The Shen EOS is based on the rela- 
tivistic mean field theory with a local density approxima- 
tion, which has been constructed to reproduce the exper- 
imental data of masses and radii o f stab le and unstable 
nuclei (see references in IShen et al.l |l998)). At the max- 
imum densities higher t han two times of sa t uration den- 
sity, muons may appear (|Wiringa et al. I ll988l ; lAkmaT et al.l 
1998). However, we neglect it, since the muon contribution 
to pre ssure at the higher density has been pointed to be very 
small (|Douchin fc Haencel|[200ll ). 

In modeling the phase transition to quark matter, there 
is a main physical uncertainty, the critical density for the 
onset of the mixed phase. Under the MIT bag model, the 
transition is determined by the value of the bag constant. 
We use the bulk Gibbs constructi on to bridge the two phases 
using the technique developed bv lGlendennind (| 19921 ). In the 
transition, two conserved quantities are the baryon number 
density (ns) and the electric charge density (: YeJiB with Y e 
being the electron fraction), 



TIB = X n B,Q + (1 - x) n B,H, 

Y e n B = X Y c,Q n B ,q + (!-%) Y c ,h n B ,H- 



(4) 
(5) 



Here \ is the volume fraction of matter in the quark phase. 
The subscripts H and Q label the number density and the 
charge fraction, Yc, in the hadronic and in the quark phase, 
respectively. Since matter in cold compact stars is in chemi- 
cal equilibrium under /3-decay with vanishing neutrino chem- 
ical potentials, the following equations have to be satisfied, 



pe Pp — pn , 

fJ-n = Uu + 2 fid, 

Up = 2u u + /id, 

Ud = Us, 



(6) 
(7) 
(8) 
(9) 



the pressure in the two phases, 



j<3 



(10) 



together with the mechanical equilibrium, the equality of 



Here for s implicity, the finite size effects on th e phase 
transition (Endo et al. 2 0061 ; iMaruvama et al.l [20061 ) are ig- 
nored. Five unknown variables to describe the mixed phase 
(x, n B,Q,nB,H, Yc,Q, Yc,h) can be solved by the five sets of 
the equations of (3,4,6,7,9). 

Figures [1] and [2] show the constructed EOS with the 
first-order phase transition. For comparison, we also plot 
the hadronic "Shen EOS" mentioned above and pure quark 
EOS("quarkB200" and "quarkB250"). The left and right 
panels are the pressure and the energy per baryon as a 
function of the baryon density, respectively. The labels of 
"mixB200(quarkB200)" and "mixB250(quarkB250)" indi- 
cate the two different bag constants, B=200 MeV fm~ 3 and 
250 MeV fm" 3 . The EOS consists of the three phases of the 
"Hadron" , "Mixed" , and the "Quark" . The critical densities 
with ni (open circles) and ni (closed circles) marking the 
boundaries between the hadron, mixed, and quark phase, 
are also given in Table [1] From the left panel and Table 1, 
it can be seen that the critical density for the mixed phase 
becomes higher for the larger bag constant. This is because 
for the large bag constant, the density should be larger to 
achieve the equilibrium between the mixed phase and the 
pure hadron phase, remembering that the pressure of the 
quark matter becomes smaller with increasing the bag con- 
stant (from Eqs. (1,2)). The higher n± for the large bag con- 
stant also leads to the higher m , the transition to the quark 
phase (Table 1). From the right panel of Figure [T] and [2] it 
is seen that the energy with the quark phase ("mixB200" 
and "mixB250") are clearly lower than the case only with 
the hadronic phase ("shen"), and also that the energy is 
generally lower for the smaller bag constant above m. 

The actual conversion process from nuclear matter to 
quark matter has been of a topic of hot de bate, such as deto- 
nation or deflagration (|Drago et alj|2007| ). The range of the 
conversion timescale is estimated to be very wide, roughly 
0.1 — 100 sec, depending upon the bag c onstant, tempera- 
ture, mass of neutron stars, and so forth (lOlesen fc Madsenl 
Il99ll ). Since our main aim of this paper is the estimation of 
the liberated energy by the conversion from neutron stars 
to hybrid stars, we do not consider the detailed combustion 
processes in this study. 

Proto-neutron stars left after supernova explosion are 
very hot (~ 50 MeV) initially, how ever, cool down be- 
low ~ IMeV in some tens of seconds (|Burrows fc Lattimerl 
1986). The newly formed neutron stars stabilize at prac- 
tically zero temperature. As mentioned, in cold neutron 
stars, the j3 equilibrium in the weak interactions can be 
well validated, with neutrinos and antineutrinos freely es- 
caping from the star. Combining the zero-temperature, zero- 
neutrino fraction, and the beta-equilibrium conditions with 
the charge neutrality condition, the thermodynamic vari- 
ables depending on the three parameters (e.g. the pressure as 
P(p, Y e ,T) with Y e being the electron fraction), can only be 
determined by a single variable (barot ropic), which we take 
to be the energy density, namely P(e) (|Shapiro fc Teukolskvl 
Il983l ) in the following. 

In the following, we introduce the baryon mass density, 
po, for convenience, defined to be the number density, n, 
multiplied a baryon mass, mo = 1.6605 x 10~ 24 g. 
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Table 1. Critical densities making the transition to 
the mixed, ni, and the quark phase, ri2, respectively. 



EOS 


ni [/fm 3 ] 


n 2 [/fm 3 ] 


mixB200 


2.30E-01 


1.29 


mixB250 


2.63E-01 


1.51 



2.2 Constructing method of equilibrium stellar 
configurations 

Employing the constructed EOS with the phase transition, 
we construct the equilibrium stellar configurations. As men- 
tioned in introduction, we pay attention to the general rela- 
tivistic and toroidally magnetized stellar configurations. The 
basic equations and t he numerical meth o ds for the purpose 
are already given in IKiuchi fc Yoshidal (|200&t ). Hence, we 
only give a brief summary for later convenience. 

Assumptions to obtain the equilibrium models are sum- 
marized as follows ; (1) Equilibrium models are stationary 
and axisymmetric. (2) The matter source is approximated 
by a perfect fluid with infinite conductivity. (3) There is no 
meridional flow of the matter. (4) The EOS for the matter 
is barotropic, which is satisfied as mentioned. (5) Magnetic 
axis and r otation axis are aligned. Because the circularity 
condition ( Wald 1983) holds under t hese assumptions, the 
metric can be w ritten in the form (|Komatsu et al I I 19891 ; 
ICook et al.lll992T ). 



ds z = -e J+p dt 2 + e za (dr z + r^dO') 
'(dip — Lodt) 2 , 



i ~f—P 2 ■ 2 

+e r sin i 



(11) 



where the metric potentials, 7, p, a, and u), are functions 
of r and 9 only. We see that the non-zero component of 
Faraday tensor F^ v in this coordinate is Fn. Integrability 
of the equation of motion of the matter requires, 



7 — Iql 



sin6»Fi2 = K(u) 



pohe 



(12) 



where K is an arbitrary function of pohe 2j r 2 sin 2 6. The vari- 
ables po and h represent the baryon density and relativistic 
specific enthalpy, respectively. Integrating the equation of 
motion of the matter, we arrive at the equation of hydro- 
static equilibrium, 



P(e) 



dP p + j 1 2^ 
— H — m( 1 — v ) 

+ P 2 2 y ' 

1 f K(u) dK , 
+ — / — ^— du = C, 
47T / u du 



(13) 



where v = (fi — uj)r sin 6e~ p with Q being the angular veloc- 
ity of the matter and C is an integration constant. Here, we 
assume the rigid rotation. The first term of equation (|13p 
depends only on EOSs. Hence, we prepare the values of the 
integral with EOS tables, precisely. This preparation enable 
us to calculate equation (|13[) precisely, though our EOSs are 
quite different from polytrope models. To compute specific 
models of the magnetized stars, we need specify the func- 
tion form of K, which determines the distribution of the 
magnetic fields. We take the following simple form, 



where b and fc are constants. Regularity of toroidal magnetic 
field on the magnetic axis requires that k ^ 1. If k > 1, 
the magnetic fields vanish at the surface of the star. In this 
study, we consider the k = 1 case because in the general rel- 
ativistic MHD simulation, IKiuchi et all (|2008cJ ) have found 
that magnetic distribution with fc 7^ 1 is unstable against 
axisymmetric perturbations. 

To solve the master equations n umerically, we em ploy 
the Cook- Shapiro- Teukolsky s cheme (|Cook et al.lll992T ) ex- 
tended by IKiuchi fc Yoshidal l|2008l ). which does not care 
about the function form of EOS. Hence, it is straight for- 
ward to update our numerical code for incorporating the 
EOS with the phase transition. 

After obtaining solutions, it is useful to compute global 
physical quantities characterizing the equilibrium configura- 
tions to clearly understand the properties of the sequences 
of the equilibrium models. In this paper, we compute the 
following quantities: the gravitational mass M, the baryon 
rest mass Mo, the total angular momentum J, the total ro- 
tational energy T, the total magnetic energy H, the mag- 
netic flux $, the gravitational energy W and the mean de- 
formation rate e, whose definitions are explicitly given in 
IKiuchi fc Yoshidal {2008). More explicitly, the mean defor- 
mation rate e is defined as 



(15) 



K(u) = bu k 



(14) 



where I xx = tt J er 4 sin 6(1 + cos 2 6)drd9 and I zz — 
2-k J er 4 sin 3 6drd9 with e being the energy density of the 
matter. Circumferential radius R c ii is defined as R c ii = 
e (7-p)/2 re Te being the coordinate radius at the stellar 
equatorial surface. 

We checked the convergence of the presented results by 
doubling the mesh numbers from the standard set of radial 
and angular direction mesh poin ts of 400 x 260. By checking 
the r elativistic virial identities (jBonazzola fc Gourgoulhonl 
Il994h for all the models, we confirm that the typical values 
are orders of magnitude 10~ 3 , and become 10~ 2 at worse 
(10~ 4 at best). These values, which is a measure of the 
numerical convergence, are almost same for the polytropic 
EOS case (Ki uchi fc Yoshidal 12008). In general, the conver- 
gence is known to become much worse for realistic EOSs, be- 
cause their density and pressure profile are not smooth due 
to phase transitions. In this respect, our numerical scheme 
works well. In Sec. 13.21 we will discuss energy releases by 
the QCD phase transition and find that the values of our 
interest are ~ 1O -2 M0, for which the numerical accuracy 
above is certificated mostly. By doubling the mesh points, 
i.e., 800(r) x 520(6*), we checked that the order of the mag- 
nitude of the released energy does not change. 

In constructing one equilibrium sequence, we have three 
parameters to choose, namely the central density p c , the 
strength of the magnetic field parameter b, and the axis ra- 
tio r p /r e . Changing these parameters, we seek solutions in 
as wide parameter range as possible to study the proper- 
ties of the equilibrium sequences. We need to calculate more 
than 50 models changing p c , b, to follow one evolutionary 
sequence for a fixed baryon mass and magnetic flux. In ad- 
dition, we change the initial angular momentum of each se- 
quence by 10 models to model the rapidly rotating to the 
non-rotating case and also 3 different EOSs of mixB200, 
mixB250, and Shen are employed. This means that we have 
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to construct at least 1500 models to explore properties of 
the magnetized compact objects with quark cores systemati- 
cally. In doing so, we use the Rosenbrock and Gram-Schmidt 
method, which is helpful to obtain the convergence of the 
solutions efficiently. 



3 NUMERICAL RESULTS 

First of all, we discuss how the EOS with the phase transi- 
tion affects the equilibrium configurations in subsection 13. II 
where we pay attention to non-rotating models. Since the 
magnetars and the high field neutron stars observed so far 
are all slow rotators, such non-rotating but highly magne- 
tized static models could well be approximated to such stars. 
Moreover the static models merit that one can see purely 
magnetic effects on the equilibrium properties because there 
is no centrifugal force and all the stellar deformation is at- 
tributed to the magnetic stress. Then in subsection 13.21 we 
move on to discuss the the releasable energy of the phase 
transition from the rotating and magnetized hadronic stars 
to the hybrid stars. 

3.1 Effect of the phase transition on the 
equilibrium configurations 

Now we discuss the equilibrium configurations of the 
non-rotating and strongly magnetized compact stars 
with/without the quark cores. Figures [3] and [4] are one ex- 
ample for the hybrid star with mixB200 EOS and for the 
neutron star with the hadronic "Shen" EOS, showing the 
distributions of the baryon density (left panel) and the mag- 
netic field (right panel) in the meridional planes, respec- 
tively. These two models have the same baryon mass of 
1.86M©, and the same magnetic flux of 2.00 x 10 30 G cm 2 , 
which mean that they are really highly magnetized models 
with the central magnetic fields of ~ 10 ls G. 

Comparing left panels of Figures [3] and [4] more con- 
centration of the matter in the center is clearly seen for the 
hybrid star (Figure[3}. The concentration is also clearly seen 
from Figure[5] Inside the inner 6 km of the core of the hybrid 
star, the mixed phase appears, leading to the enhancement 
in the compression of the matter due to the softening of the 
EOS (e.g., Figures 1 and 2). 

From right panel of Figure [4] the strong toroidal field 
lines are seen to behave like a rubber belt, wrapping around 
the waist of the neutron star. It is found that the magnetic 
fields frozen-in to the matter, are also compressed by the 
presence of the quark phase for the hybrid star (right panel 
of Figure [3}- In fact, the pinching of the field lines in the 
panel corresponds to the surface of the quark core, ~ 6km 
in radius. These qualitative features are also true for the 
other equilibrium models. 

In addition, a general trend in the equilibrium configu- 
rations for the hybrid stars is that they are more compact 
than the neutron star, due to the softness of the EOS. Given 
the same stellar baryon mass (Mo), the gravitational mass is 
smaller up to ~ O.OIA/© than for the neutron stars, reflect- 
ing their smaller energy of the quark matter (right panel of 
Figure [TJ . 

Given the fixed magnetic flux, we can construct equi- 
librium sequences by changing the central density. To char- 



acterize the features of the hybrid stars, we pay attention to 
the model with the maximum mass (dM /dpo.c = 0) along 
each sequence, which we call as the maximum mass sequence 
for convenience. In table 2, the important physical quanti- 
ties are summarized for the maximum mass sequences with 
different magnetic fluxes. It is here noted that the maximum 
value of the magnetic flux in the table ($ = 2.5 x 10 30 G cm 2 ) 
corresponds to the non-convergence limit, beyond which any 
solutions cannot converg e with the present numerical scheme 
jKiuchi fc Yoshida| [2008). Albeit with such limitations, the 
field strength is already enough high to affect the configura- 
tions and we can well study the magnetic effects on them. 

Since the compression of the matter is more enhanced 
for the smaller bag constant, the hybrid stars with the 
smaller bag constant become more compact (see Rdr Ta- 
ble 2). The maximum magnetic fields are found to become 
up to about 30 % larger for the hybrid stars than for the 
neutron stars (compare -B max for mixB200 with for Shen). 
It is also found that the hybrid stars with the smaller bag 
constant become more prolate (smaller values of e) and also 
their maximum masses become smaller up to ~ 10% than 
those for the larger bag constant models. All these features 
are helpful to understand the properties of the evolution 
tracks of the hybrid stars, which we discuss from the next 
section. 

3.2 An evolutionary track to a hybrid star 

Based on the equilibrium configurations mentioned above, 
we now move on to discuss an evolutionary path to the for- 
mation of a hybrid star due to the spin-down of magnetized 
and rotating neutron stars. 

The evolution scenario we have in mind is illustrated 
in Figure [5] Let us consider evolution of a single pro- 
toneutron star, left after core-collapse supernova explosion, 
which rotates with the mass-shedding limit ((A) in Fig- 
ure [5} . The maximum magnetic fields deep inside the core 
are taken to be ~ 10 18 G, which could be sustained due 
to a — Q dynamos in such a rapi dly rotating neutron star 
i|Thompson fc Duncar]|l993l . I1996T I. It should be noted that 
the possibility of such ultra-magnetic fields has not been 
rejected so far, because what we can learn from the obser- 
vations of magnetars by their periods and spin-down rates 
is only their surface fields (~ 10 15 G). During their evolu- 
tion, the baryon mass and the magnetic field strength are 
assumed to be constant for simplicity. Such models may 
model the evolution of the isolated compact stars, losing 
angular momentum via the gravitational radiation and/or 
magnetic breaking (from (A) to (B)). The phase transition 
to the hybrid stars is expected to take place during the evo- 
lution (shown as (B) to (C) in Figure|SJ). At the moment, the 
baryon rest masses and the angular momenta are assumed to 
be conserved, which may be justified because the timescale 
of the conversion, albeit uncertain of 0.1-100 s, are too short 
compared to typical evo lutionary timescale of co mpact stars 
(more than 10 years) (01 cscn fc Madsenlll99ll ). After the 
transition, the newly born hybrid star evolves, again losing 
the angular momentum, to settle down to the magnetized 
and no-rotating hybrid star finally (see (C) to (D)). 

Now using the left panel of Figure [7] we proceed to dis- 
cuss quantitatively the evolution tracks mentioned above. 
The red solid the baryon mass 'Mo = I.86M0' rotating 
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Table 2. Global physical quantities for the maximum gravitational mass models of the constant magnetic 
flux sequences of the non-rotating stars. The order of EOSs from top to down indicates for their hardness 
from the hadronic "Shen EOS" to the softest "mixB200" EOS. 



$[10 30 G 


cm 2 ] 


P0,c[10 15 g/cm 3 ] 


M[M Q ] 


M [M Q ] 


-Rcir[km] 


Bmax[10 


18 G] 


H/\W\ 


e 


Shen 


0.00E-4 


-00 


1.16E+00 


2.46EH 


-00 


2.86E+00 


1.35E+01 


0.00E4 


-00 


0.00E4-00 


0.00E4-00 


1.50E4 


-00 


1.01E+00 


2.45E^ 


-00 


2.83E+00 


1.39E+01 


7.34E-01 


3.81E-03 


-7.40E-03 


2.00E4 


-00 


1.05E+00 


2.45E4 


-00 


2.84E+00 


1.39E+01 


9.83E-01 


6.59E-03 


-1.30E-02 


2.50E4 


-00 


1.17E+00 


2.46EH 


-00 


2.85E+00 


1.36E+01 


1.26E4 


-00 


9.48E-03 


-1.86E-02 


mixB250 


0.00E4 


-00 


1.18E+00 


1.83Eh 


-00 


2.03E+00 


1.47E+01 


O.OOE-t 


-00 


0.00E4-00 


0.00E4-00 


1.50E4 


-00 


1.17E+00 


1.83E4 


-00 


2.02E+00 


1.40E+01 


9.38E-01 


8.28E-03 


-1.97E-02 


2.00E4 


-00 


1.18E+00 


1.83E4 


-00 


2.02E+00 


1.41E+01 


1.23E4 


-00 


1.45E-02 


-3.50E-02 


2.50E4 


-00 


1.64E+00 


1.82E4 


-00 


2.00E+00 


1.33E+01 


1.76E4 


-00 


2.02E-02 


-4.63E-02 


mixB200 


0.00E-4 


-00 


1.44E+00 


1.70E4 


-00 


1.88E+00 


1.30E+01 


0.00E4 


-00 


0.00E4-00 


0.00E4-00 


1.50E4 


-00 


1.29E+00 


1.69E4 


-00 


1.86E+00 


1.34E+01 


1.08E4 


-00 


9.62E-03 


-2.25E-02 


2.00E4 


-00 


1.35E+00 


1.69E^ 


-00 


1.86E+00 


1.34E+01 


1.44E4 


-00 


1.65E-02 


-3.92E-02 


2.50E4 


-00 


1.65E+00 


1.69E^ 


-00 


1.86E+00 


1.31E+01 


1.86E4 


-00 


2.39E-02 


-5.59E-02 



with the mass shedding limit (corresponding to (A) in Fig- 
ure [6]). It is noted that the choice of the baryon mass 
(Mo = I.86M0) is determined by the maximum mass of the 
non-rotating hybrid star with the mixB200 EOS as discussed 
in subsection 13.11 and also that this baryon mass constant 
sequence has the largest releasable energy at the transition 
among all the models, as we show this point later. There are 
many possible paths to convert to hybrid stars with differ- 
ent masses. However, we focused on the maximum energy 
release with QCD phase transition in this paper. It is pos- 
sible to select a single equilibrium sequence by keeping the 
baryon rest mass and the magnetic flux constant simulta- 
neously. Thus the evolution of the protoneutron star can 
be described along the line (pink dotted line), going right 
from one point on the red solid curve. As mentioned, our 
models assume that all equilibrium sequences begin at their 
mass-shedding limits and continue to non-rotating equilib- 
rium hybrid stars, at which the sequences end. The final 
fate is shown by the green dashed line in Figure [7] which is 
the non-rotating hybrid star with the same baryon mass of 
'M = 1.86M ', here for the mixB200 EOS (corresponding 
(D) in Figure!}. 

Between pi and pi (the vertical lines in Figure [7]), the 
phase transition should take place, producing the jump to 
the sequence of hybrid stars (from (B) to (C) in Figure |6]). 
To see clearly, the transition is shown as the vectors from cir- 
cle (o) to filled circle(») in the right panel of FigureUJ which 
is just the magnification of the left panel near the phase tran- 
sition. Here, 'J = 1.00' ('J = 0.70', 'J = 0.50', 'J = 0.25') 
denotes the angular momentum in unit of 10 49 g cm 2 s _1 . We 
estimate the energy release at the phase transition simply 
by the difference in the mass-energies between the neutron 
star and the hybrid star, 

AAf = M NS (M , J) - M HS (M , $, J), (16) 

with conserving $ and J through the transition. 

In Figure [8] we show the releasable energy of AM as a 



function of the angular momentum for the mixB200 EOS. 
Top panel is for the baryon mass of Mo = I.86M0, which 
is a special evolutionary track in the sense that the baryon 
mass is the maximum among non-rotating models, so that 
the releasable energy is largest among the computed models. 
Also for the other baryon mass models (middle and bottom 
panels), it is seen that the energy release becomes smaller 
for the larger angular momentum. This is simply because 
stronger centrifugal forces prevent the compression (see the 
central density (p c ) in Tables [3] and [4} , making the quark 
cores smaller, which weakens the effect of the phase transi- 
tion. In all the computed models, the maximum energy re- 
lease is found to be < 0.01A/© . This energy is equivalent to < 
2 x 10 52 erg, which is truly larger compared with the energies 
of supernovae (~ 10 51 erg), however much smaller than those 
previously estimated (~ 10 53 erg) i n which the conversion to 
the strange quark star is assumed rtGourgoulhon et al 1 ll999l : 
I Yasutake et al.ll2005l ; IZdunik et al.ll2007t ). We note that AM 
becomes less than AM < 1O _3 M0 for the larger bag con- 
stant models of mixB250. This is because the large bag con- 
stant suppresses the conversion from hadron matter to quark 
matter as already mentioned. Unfortunately, exact values 
of AM for the mixB250 models cannot be found by the 
present scheme because AM becomes more than 3 orders- 
of-magnitude smaller than M, while M can be estimated 
at most with the numerical accuracy of 10~ 2 for this case. 
This is also the case for small baryon mass cases (see the bot- 
tom panel of Figure 8). For getting more precise solutions, 
we may need to employ th e so-called spectral me thod. Al- 
though the LORENE code l|Bonazzola et al.lll998t ) employs 
the method, it cannot treat the magnetic fields yet. The 
implementation of the method is still a major undertaking, 
which we pose as a next task of this paper. 

Other numerical values characterizing the phase tran- 
sition are given in Tables [3] and [4] Table [4] is for the non- 
magnetized models, which is given for comparison. For mod- 
els with smaller angular momentum (J in the tables), it can 
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be readily seen that the gravitational masses (M) and the 
circumferential radii (-R c ir) become smaller, while the cen- 
tral baryon densities (po.c) and the maximum strength of 
the magnetic fields (£> ma x) become larger, all simply due 
to the smaller centrifugal forces. If the newly-born neutron 
star evolve to the non-rotating neutron star without experi- 
encing the conversion, the increase of the central density is 
less than ~ 20% (compare po,c of " Shen" at the mass shed- 
ding limit with the one at no-rotating limit in Table [3J. On 
the other hand, the maximum density is found to increase 
about several times larger by the phase transition (compare 
po, c between "Shen" with "mixB200" in Table |3J|. Frozen-in 
to the matter, the maximum magnetic field increases up to 
50% while the change is an order of unit percent in absence 
of the conversion (see Figure O. Due to the conservation 
of angular momentum, the angular velocity increases about 
Afi/fi ~ 30% after the conversion (see Figure PTOl and Ta- 
bles [3). Also in this case, it is noted that the spin-up is sup- 
pressed for the larger bag constant models (see Table 3). We 
speculate that such spin-up, much larg er than the canonical 
fitch es of magnetars (Afi/fi ~ 10~ 7 . [Woods fc Thompson] 
2006)), could be new observable signatures, marking the 
occurrence of the phase transition. 



4 SUMMARY AND DISCUSSION 

In this study, we investigated structures of general rela- 
tivistic hybrid stars containing super-strong magnetic fields. 
Pushed by the results of recent stellar evolution calcula- 
tions and the outcomes of recent MHD simulations of core- 
collapse supernovae, we treated the toroidal fields only. Us- 
ing the bulk Gibbs construction, we modeled the EOS with a 
first order transition by bridging the MIT bag model for the 
description of quark matter and the nuclear EOS by Shen 
et al. We found that the presence of the quark phase can af- 
fect the distribution of the magnetic fields inside the hybrid 
stars, leading to the enhancement of the field strength about 
30 % than for the normal neutron stars. Using the equilib- 
rium configurations, we explored the possible evolutionary 
paths to the formation of hybrid stars due to the spin-down 
of magnetized and rotating neutron stars. For simplicity, the 
total baryon mass and magnetic flux are taken to be con- 
served during the evolution but that the angular momenta 
are lost gravitational waves and/or magnetic breaking. We 
found that the energy release by the conversion to hybrid 
stars is typically < 10 52 erg, smaller than previously esti- 
mated for the conversion to the strange quark stars. 

In association with the vast energy release, it is natural 
to expect the emissions of gravitational waves. Amplitude of 
the gravitational wave s associated w ith the phase transition, 
h, is given as follows (|Pachecolll998h . 



damping tim escale and a is ~ 100 ms an d 10 4 , respectively 
inferred from lAbdikamalov et al.l (|2008f l. The resulting am- 



1 

ruj 



(17) 



Here, -Etot, a, T, r, and ui are energy release, fraction of the 
energy release emitted in form of the gravitational waves, 
damping time scale of the gravitational wave, distance to 
the source, and angular velocity, respectively. To set an ab- 
solute upper bound, we could choose E to t ~ 10 52 erg and 
(jj ~ 100 rad s _1 according to Table [3] We assume that the 
source is in our galactic center of r ~ 10 kpc and that the 



plitudes become as high as h ~ 10 



peaking at ~ kHz, 



which are surely within the detection limits for the laser in- 
terferometers on line such as LIGO, VIRGO, GEO600, and 
TAMA300. Considering the optimal sensitivity of the in- 



terferometers of h 



10" 



at ~ 1 KHz, such signals, far 



stronger th an the ones from can oni cal core-collapse supe r- 
novae (e.g.. iKotake et all (|2009al lbh: iMurphv et all (j2009h ^. 
are possibly visible for a source out to Megaparsec distances. 
In combination with the signatures of the spin-up discussed 
before, we speculate that the waveforms of such strong 
gravitational waves could provide us hints of properties 
of QCD phase transition. To find the accurate waveforms, 
general relativistic simulations are indispensable (Kiuchi et 
al. in preparatio n). Moreover accura t e predictions for neu- 
trino emissions l|Gentile et alj 1 19931 ; iNakazato et alj 120081 ; 
ISaeert et al J l200q) and their de t ectab ility are remained to 
be studied (e.g., Kawagoe et al.l (|2009f )). This paper should 
help to construct the initial conditions for such studies. 

lAlcock et al.l (|l986l ) claimed that the conversion from 
the hadronic to the mixed phase occurs instantaneously in 
the weak interaction timescales (such as d + u «-> s + u) 
when the central density exceeds the critical value of pi in 
Figures 1 and 2. If so, the newly-born neutron stars may 
convert to the hybrid stars in the postbounce phase, be- 
fore they evolve from A to B illustrated in Figure [6] In 
this case, the released energy may be used to power orig- 
inal core-collapse supernovae or more energetic supernovae 
such as hypernovae. Although the formation paths to hybrid 
stars are different from the ones discussed in this paper, the 
estimation method of the energy release here are also ap- 
plicable if we implem ent the EOS at finite temp erature and 
high lepton fraction (jYasutake fc K ashiwa 2009), which we 
plan to study as a sequel to this paper. 

It has been pointed out that the super-strong magnetic 
field more than 10 18 G can affect the stiffness of the EOS 
l|Broderick et al.ll200ol ). As mentioned, the maximum mag- 
netic field increases about ~ 30% by the phase transition, 
while it increases only unit percent if the compact stars 
evolve to the non-rotating ones in absence of the transition. 
This suggests that the magnetic effects should be seriously 
taken into account to our models t reating the phase transi- 
tion. E s pecially for the quark phase. iFukushima fc Warringal 
i|2008l ); lNoronha et all l^OOTtl have recently found that the 
energy gaps of magnetic color-flavor-locked phase are oscil- 
lating functions of the magnetic field. Such effects are re- 
mained to be studied. 

We should comment about the very high values of the 
magnetic fields, larger than 10 18 G. iRuderman et al.l (120001 ) 
suggested that such high magnetic fields can be brought to 
the stellar surface by buoyancy forces. According to their 
result, we estimate the buoyancy time scale in our models. 
The buoyancy force can be expressed as 



87TC 2 



(18) 



where c s and g are the sound speed and the gravitational 
acceleration, respectively. Putting typical values inside com- 
pact objects into the equation, we roughly estimate the time 
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Table 3. Global physical quantities for the equilibrium sequences of the rotating stars with Mq = 1.86Mq, $30 = 2.00. 



EOS 


p ,c[10 15 g/cm 3 ] 


M[M Q ] 


R cir [k: 


Hi] 


n [10 3 rad/s] 


T/\W\ 


H/\W\ 


e 


Bmax[10 18 G] 


Shen(MS) 


4.11E-01 


1.73E+00 


2.22E4 


-01 


1.59E+00 


9.41E-02 


2.89E-02 


2.00E-01 


9.61E-01 






J = 1.00 


X 10 49 


g cm 












Shen 
mixB200 


4.51E-01 
5 11E-01 
8.26E-01 


1.72E+00 
1.71E+00 


1.70E4 
1.61E4 


-01 
m 
-01 


8.73E-01 

Q Q1 TT HI 

1.30E+00 


3.69E-02 
3.78E-02 


2.47E-02 

9 49F 09 
2.11E-02 


5.40E-02 

K AOF 09 

7.01E-02 


9.65E-01 

Q fiSTT 01 

1.10E+00 






J = 0.75 


x 10 49 


g cm 












Shen 

mivR9 c; in 

mixB200 


4.63E-01 
9.63E-01 


1.71+00 
1.71E+00 


1.64E4 
1.51E4 


-01 
01 
-01 


6.69E-01 

7 AOTT 01 

1.10E+00 


2.15E-02 

9 1 09 
2.24E-02 


2.37E-02 
9 ^'^F 09 

1.93E-02 


7.26E-03 

O.OZIZj-UO 

2.74E-02 


9.70E-01 

Q 01 

1.20E+00 






J = 0.50 


x 10 49 


g cm 












Shen 
mixB250 
mixB200 


4.72E-01 
5.93E-01 
1.08E+00 


1.71E+00 
1.71E+00 
1.70E+00 


1.60E4 
1.59E4 
1.43E4 


-01 
-01 
-01 


4.54E-01 
5.16E-01 
8.07E-01 


9.83E-03 
9.77E-03 
1.04E-02 


2.31E-02 
2.24E-02 
1.81E-02 


-3.05E-02 
-2.82E-02 
-8.49E-03 


9.70E-01 
9.75E-01 
1.31E+00 






J = 0.25 


x 10 49 


g cm 












Shen 
mixB250 
mixB200 


4.79E-01 
6.15E-01 
1.21E+00 


1.71E+00 
1.71E+00 
1.69E+00 


1.58E4 
1.57E4 
1.42E4 


-01 
-01 
-01 


2.28E-01 
2.60E-01 
4.39E-01 


4.53E-03 
4.43E-03 
2.82E-03 


2.28E-02 
2.20E-02 
1.71E-02 


-4.82E-02 
-4.57E-02 
-3.15E-02 


9.70E-01 
9.75E-01 
1.39E+00 


no rotation 


Shen 
mixB250 
mixB200 


4.86E-01 
6.36E-01 
1.35E+00 


1.71E+00 
1.71E+00 
1.69E+00 


1.57E4 
1.55E4 
1.34E4 


-01 
-01 
-01 


0.00E+00 
0.00E+00 
0.00E+00 


0.00E+00 
0.00E+00 
0.00E+00 


2.25E-02 
2.16E-02 
1.65E-02 


-6.35E-02 
-6.03E-02 
-3.92E-02 


9.72E-01 
9.75E-01 
1.44E+00 



Table 4. Same as Table |3] but for * 30 = 0. 



EOS 


P0,c[10 15 g/cm 3 ] 


M[M Q ] 


R c i T [km] 


Q [10 3 rad/s] 


T/\W\ 


e 


Shen(MS) 


3.81E-01 


1.74E+00 


2.10E+01 


1.45E+00 


1.12E-01 


2.79E-01 






J = 1.00 


x 10 49 g cm 


2 /s 






Shen 
mixB250 
mixB200 


4.40E-01 
4.62E-01 
7.80E-01 


1.72E+00 
1.72E+00 
1.71E+00 


1.66E+01 
1.66E+01 
1.58E+01 


8.08E-01 
8.30E-01 
1.18E+00 


3.58E-02 
3.54E-02 
3.72E-02 


1.07E-01 
1.06E-01 
1.11E-01 






J = 0.75 


x 10 49 g cm 








Shen 
mixB250 
mixB200 


4.51E-01 
5.16E-01 
9.02E-01 


1.71E+00 
1.71E+00 
1.70E+00 


1.61E+01 
1.61E+01 
1.49E+01 


6.23E-01 
6.69E-01 
1.01E+00 


2.09E-02 
2.07E-02 
2.18E-02 


6.49E-02 
6.43E-02 
6.90E-02 






J = 0.50 


x 10 49 g cm 


2 /s 






Shen 
mixB250 
mixB200 


4.61E-01 
5.51E-01 
1.03E+00 


1.71E+00 
1.71E+00 
1.70E+00 


1.58E+01 
1.57E+01 
1.42E+01 


4.25E-01 
4.69E-01 
7.51E-01 


9.54E-03 
9.49E-03 
1.02E-02 


3.06E-02 
3.04E-02 
3.33E-02 






J = 0.25 


x 10 49 g cm 


2 /s 






Shen 
mixB250 
mixB200 


4.67E-01 
5.74E-01 
1.14E+00 


1.71E+00 
1.71E+00 
1.69E+00 


1.56E+01 
1.55E+01 
1.37E+01 


2.15E-01 
2.42E-01 
4.08E-01 


2.54E-03 
2.52E-03 
2.65E-03 


8.48E-03 
8.46E-03 
9.10E-03 



no rotation 



Shen 4.75E-01 1.71E+00 1.55E+01 0.00E+00 0.00E+00 0.00E+00 

mixB250 5.82E-01 1.71E+00 1.54E+01 0.00E+00 0.00E+00 0.00E+00 

mixB200 1.19E+00 1.69E+00 1.35E+01 0.00E+00 0.00E+00 0.00E+00 
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scale t as 



10" 3 - 10" 



1/2 f_^y /2 f-^—\ f^-y 1 ( 9 v 1/2 (i 9 ) 

10 15 g/cm 3 I UO B cm/ I 10 10 cm/s J\ 10 18 G/ I 10 14 cm/s 2 / 



where R is the typical size of compact stars. Clearly such 
magnetic fields are not stable against the buoyancy forces 
in stellar evolution time scale ~ 10 3 y r. Based on the MHD 
simulations in full general relativity, iKiuchi et al.l (|2008c) 
recently showed that the toroidal configurations of neutron 
stars is really dynamically unstable due to the buoyant in- 
stability. But more important finding relevant to this study, 
is that the toroidal magnetic fields settle down to a new 
equilibrium state with the circular motions in the meridian 
plane. In the new equilibrium configurations, the toroidal 
fields almost equivalent to the strength before the onset of 
the instability, are expected to be much stronger that the 
poloidal ones. This suggests that our models presented here, 
albeit unstable to the buoyant instability, could be helpful 
to understand equilibrium configurations of magnetized hy- 
brid stars. Such stability analysis is an important issue yet 
to be studied. 

Our treatments of the phase transition should be more 
sophi sticated. We plan to employ the so-called NJL model 
(e.g., lYasutake fc Kashiwal (J2009)) instead of the simple 
MIT bag model. Moreover we sha l l tak e into account fi- 
nite size effects l|Maruvama et af1 I2006T ) such as quark- 
nuclear pasta structures. We considered cold compact ob- 
jects, namely with zero-temperature and zero-neutrino frac- 
tion. We plan to extend this study to the finite temperature 
with the non-zero neutrino fraction, which should be useful 
for studying earlier evolutions of compact stars soon after 
their formation. 

The evolution sequence we explored is a very simplified 
one. For more realistic estimations, one needs the two di- 
mensional evolutionary calculation in full general relativistic 
framework, not an easy job, due to the treatment of convec- 
tion, combustion, cooling and/or heating processes, nucle- 
osynthesis on the surf ace, and so f o rth. A pplying the method 
recently reported by iPons et all i|200Sl ). we think that we 
could study the magneto-thermal evolution of neutron stars 
to the hybrid stars. By doing so, we hope to investigate 
the peculiar properties in the light curves, whic h has been 
pointed out to be another observationally visible ([Page et al.l 
l200d ; iBlaschke et aHl200Cl . l200ll ; iGrigorian et alj|2005h . All 



of such studies should be indispensable to pave the way 
for the understanding of the long-veiled phase-transition 
physics from the astrophysical phenomena. 
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Figure 1. The constructed EOS under bulk Gibbs condition. The labels of "quarkB200" and "mixB200" indicate a pure quark matter 
EOS and an EOS with mixed phase with B = 200 MeV fm -3 , whereas "Shen" indicates Shen EOS of pure hadronic matter. Left and 
right panels are the pressure and the energy per baryon as a function of the baryon number density, respectively. 
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Figure 2. Same as Figure 1 but for B = 250 MeV fm 
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Figure 3. Distribution of (a): baryon density [g/cm 3 ] and (b): magnetic field with "mixB200" EOS. Here, the central baryon density 
PO.max = 1-35 X 10 15 g cm -3 , the baryon mass Mo = 1.86A/q, the gravitational mass M = 1.69M0, the flux normalized by units of 
10 30 G cm 2 , <I>3o = 2.00, the ratio of the magnetic energy to the gravitational energy, ff/|W^| = 1.65 X 10 -2 , maximum magnetic fields, 
Bmax = 1.44 X 10 18 G. 
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Figure 4. Same as Figure l4l but for hadronic "Shen" EOS, Mq = 1.86Mq star with the central baryon density po,max = 4.86 X 10 14 
cm" 3 , M = 1.71M , H/\W\ = 2.25 x lCT 2 , B m ax = 0.97 x 10 18 G. 
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Figure 5. Distribution of baryon density for Figure \3\ and Figure [4] 
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Figure 6. Schematic drawing, showing our scenario of conversion from neutron stars (NSs) to hybrid stars (HSs). The baryon rest mass 
and the magnetic flux are assumed to conserve for modeling the isolated neutron stars that are adiabatically losing angular momentum 
via the gravitational radiation and/or magnetic breaking. We furthermore assume the conservation of the angular momentum at the 
conversion. 
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Figure 7. Gravitational masses versus central effective baryon densities at the constant magnetic flux, <I>3o = 2.00. The right panel is 
magnified figure of the left panel. The solid line of MS(Shen) is for mass shedding limits of neutron stars. The dashed line of NR(mixB200) 
is for no-rotating hybrid stars with "mixB200" EOS. The line of 'Mo = 1.86 (Shen)' is one example in all baryon mass constant lines 
for neutron stars. 'Mo = 1.86 (mixB200)' is same as 'Mo = 1.86 (Shen)', but for hybrid stars. Conversions from neutron stars to hybrid 
stars are shown as the vectors from circle (o) to filled circle (•). Here, J is the angular momentum in unit of 10 49 g cm 2 s — 1 . The lines 
of pi and p2 are critical densities of phase equilibrium for hadron phase and quark phase shown in Figure 1. and Table 1. 
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= 2.00) and non-magnetized models ($30 = 0) for Mo = 1.86Mq (the upper panel), Mq = 1.81Mq (the middle panel), and 
Mo = 1.76Mq (the lower panel). Note here that these are the case for the mixB200 EOS. For the mixB250 EOS, exact values of AM 
cannot be estimated due to the limitation of the present scheme. Unfortunately it is also seen here, for example, for the case of smaller 
baryon mass for J > 0.4 X 10 49 [g cm 2 s — x ] (bottom)(see text for more details). 
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Figure 9. Central baryon density versus maximum magnetic field for same stars of Figure [7] 




Figure 10. Central baryon density versus angular velocities for same stars of Figure [7] Angular velocities are monotonically decreased 
without the conversion, shown as '>>'. 
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